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Abstract 

Dynamics of a fiuxon in a stack of coupled long Josephson junctions is studied 
numerically. Based on the numerical simulations, we show that the dependence of the 
propagation velocity c on the external bias current 7 is determined by the ratio of the 
critical currents of the two junctions J. 

1 Introduction 

The Stacked-Josephson junctions we consider here consist of three slabs of superconducting 
material which have an insulating barrier between two superconductors. An important 
application of these long Josephson junctions (LJJ) is the magnetic flux quanta (fluxons). 
The fluxon, which is a circulating current, can move along the junctions if there is a biased 
current applied. 

The one dimensional stacked-LJJ can be modelled by a perturbed sine-Gordon equation 
which, in normalized form, may be written as 

<l>lx - 4>u - sin 0^ - S4>1^ = a4)\ - 7, , , 

- 4>l - sin 4>VJ - = - 7- 

Here is the superconducting phase difference across the junction i pp. The term a4>t 
models the damping due to quasi particle tunnelling, a > 0. The driving force, corresponding 
to the normalized bias current 7, typically lies in the interval < 7 < 1. S' G [— 1,0] is the 
coupling parameter of the two junctions and J is the ratio of the critical current of of the 
two junctions. 

One of the solutions of the unperturbed-sine-Gordon equation {(f>xx — 4'tt = sin 0) is the 
soliton 

00 = 4arctan [exp ( — ° )], (2) 

Vl — c"^ 

which is moving with velocity c, |c| < 1. In the unperturbed equation this parameter is not 
determined. The wave speed will become a function of 7 when a and S do not vanish and 
the damping is balanced by the driving force. Therefore, c is a function of 7 for given a 
and S. The free parameter ^0 determines the initial position of the fluxon. In physics, this 
soliton solution is related with units of the flux quantum or fluxon. The system has, 
when S* = a = 7 = 0, the solution (0O: 0) in the primary state [1|0]. 

In this paper, we investigate the relation of c and 7 for spesific values of a and S in 
the primary state [1|0]. To represent the configuration, the notation [n\m] is used to denote 
the state where n solitons are located in the first junction and m solitons in the other. 
Negative numbers are used for antisolitons. In Sec. [21 we consider the case of identical 
junctions. Unequal junctions are considered in Sees. |21 and 0] The junctions considered in 
these sections differ only on the critical currents. In Sec. [SI we draw conclusions. 
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2 Identical Junctions 



We look at travelling wave solutions to the equation So we assume that the solution 
only depends on £^ = x — ct. The partial differential equation is thus reduced to a system of 
ordinary differential equations 

(1 - c2)0i - sin 0i - S(t)l^ = -ca^ - 7, 

(1 - c2)02^ - sin 02/ J - S(l)l^ = -ca(/.| - 7. ^ 



0.6- 



0.4- 




Figure 1: The relation between 7 (vertical axis) and c (horizontal axis) which is known 
as the IV-characteristic is presented. Picture (a) shows the complete relation of 7 and c. 
Picture (b) magnifies a part of the plot (a) where it becomes clear that the curve indeed 
contains a spiral. 

Considering ^ in {(p^, (/)|, 0|)"^ phase space, the problem of computing solitary wave 
profiles and speeds for with [1|0] state is equivalent to finding heteroclinic orbits con- 
necting two fixed points: 

(0\02) = (arcsin7,O) (4) 
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and 



, 0^) = (arcsin7 + 27r, 0). 



(5) 



We can calculate the parameter combinations for which there exist such a heteroclinic solu- 
tion numerically, using the program AUTO First, we consider the case when J = 1. In 
this case, the system of equations (jSJ can be rewritten to 



S 



s 
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sm( 



■ sm ( 



r(ac(/)?'^ - 7) + ac(/)P - 7 (6) 



from which we see a singularity at 1 



= ±S or |c±| — \J\ ± S. These values of c are 



called the Swihart velocity. Note that |c+| < |c_| because S is negative. 

In Fig. n the numerically obtained (c, 7) picture (in the future we call these kind of 
pictures IV-characteristics) is shown. The numerical values used in all calculations are 
S = -0.25 and a = 0.17. 

Somehow, near the edge of the IV-characteristic, we have a spiral shaped curve. Previ- 
ously, this kind of spiral in the IV-characteristic has been observed in a single long Josephson 
junction by Brown et al. 4 . Later on, it is shown in [101 that such a spiral exists in 
the single Josephson junction if the dissipation due to surface resistance, i.e. adding a term 



w^^j in Eq. is taken into account. 



However, there is a difference between the spiral we have in this case and spiral in a 
single long Josephson junction. The spiral in the single long Josephson junction has its 
center located at |c| = 1 which is the maksimum velocity possibly achieved by a fluxon. In 
our case here, the center is below the lowest Swihart velocity |c+| = ^Jl -f S 'i . The first 
four turning points are given in Tabled The plot of the solitons related to those turning 
points is given in Figure |2| 
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c 


0.7489735 


-0.8622143 


0.4900705 


-0.8625241 


0.5161983 


-0.8625433 


0.5121004 


-0.8625403 



Table 1: The first four turning points of the curve in the 7 — c space. 

Goldobin et al. |S] give a 'proof which shows that the fluxon in the state [1|0] with 
J = 1 cannot pass the lowest Swihart velocity |c+|. The result says that if jc+l is reached, 
there is at least a fluxon in the second junction which then gives the [1|1] state. 

3 Unequal Junctions: J > 1 

An interesting situation occurs if the two junctions forming the stack are not identical. 
We assume that only the critical currents of the two are unequal, which translates to the 
condition that J in Eq. (O is not equal to 1. 

Goldobin et al. |H] mention for the first time that there is a 'back-bending' phenomenon 
in the IV-characteristic. With AUTO we can continue the backbending such that we get 
the full-branch of the IV-characteristic. The result is shown in Fig. 13 

Goldobin et al. |S] also show that the velocity jc+j cannot be reached at this case. 

In Fig. 01 two solutions for the same value of 7 are plotted. In the first junction there 
is hardly any change, but in the second junction there is a big difference. On the left 
branch there is only a small 'image' in the second junction, which then if we follow the 
IV-characteristic, this image has grown to a large hump. 
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The stability of solutions on the right branch has been tested numerically in (TT|. By 
taking the solution as calculated with AUTO and using this as initial solution to the PDE- 
solver we conclude that all solutions on the right branch are unstable. All solutions at 
the left branch up to the point where drjjdc — are stable. Solutions on the right branch 
are in the domain of attraction of the solution for the same value of 7 on the left branch. 

The instability of the solution can be seen clearer when we follow the right branch in 
the direction of decreasing 7. In the limit 7 tends to 0, there is a fluxon-antifluxon pair in 
junction 2 or the stack is in — 1] jS]- Because the fluxon and antifluxon are attracting 
each other, if we apply a bias current, they will collide and a new different state is formed. 
This is the physical argument to show that this state is unstable. 

4 Unequal Junctions: J < 1 

A somewhat different behavior happens when we look at case the J < 1. While in the case 
J > 1, the velocity is bounded by c-|_, |c| < |c+|, in the case J < 1 we have that |c| can 
exceed the value |c+|. When fluxons move with velocities above the lowest Swihart velocity, 
Cherenkov radiation takes place. 

Cherenkov radiation is the emission of waves behind the moving fluxon. Behind the 
fluxon we will find an oscillating tail which is the emitted waves. E. Goldobin et al. [HI [7| 
observe Cherenkov radiation happening in the [1|0] state numerically and experimentally. 

The IV-characteristic of the state is given in Fig.l^a). Unfortunately, the IV-characteristic 
stops at c+ because our equation becomes singular at that value and AUTO cannot con- 
tinue the calculation through the singularity. Therefore, all solutions that lie in this IV- 
characteristic do not have such emitted wave behind the fluxon. 

If we keep increasing the applied bias current, the velocity of the fluxon increases up to 
the maximum velocity |wi\ai;|. This maximum velocity is a function of J, S, a By now, 
the analytic form of the maximum velocity for this state is still unknown. 

Another branch is found when we use one solution on the right branch of the IV- 
characteristic for J > 1 as the initial solution of AUTO. Using this initial condition, then 
we decrease the parameter J up to J = 0.5. In fact, we get a backbcnding curve for the 
IV-characteristic of the case J < 1. The curve is given in Fig. El^b). If for J > 1, we see 
hardly any changes when comparing solutions in the first junction, while something differ- 
ent happens in the second junction (Fig. QJ, for the case we have the opposite result. In 
junction 2 we almost have seen nothing, while there is a difference in the solutions in the 
first junction. In Fig. |Blwe compare two solutions on the left and right branch for the same 
value of 7. 

A further analysis is needed to see the stability of solutions lying on this branch. We 
conjecture that all the solutions are unstable. The same argument might be used as in 
the case J > 1 before. Numerical simulation confirms as well that the solutions go to the 
corresponding solutions at the other branch. We have done that also using PDE-solver 
where we use a solution resulted by AUTO as the initial solution of the PDE-solver. 

5 Conclusion 

We have given several calculations showing that stacked-Josephson junctions have very rich 
behavior. Just with [1|0] state and changing J, we already have qualitatively different 
IV-characteristics, while there are many more states that can exist in the junctions. 

How the IV-characteristics corresponding to the state [n\m] change when J passes 
through 1 needs to be investigated further. Also the stability of the various patterns needs 
to be determined. 
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Goldobin et al. [HI state that for the [1|0] state not only the condition J < 1 leads to 
|c| > |c+|, but also the asymmetry of the appUed bias ern'rent of the two-fold stack, i.e. 
when we have 7^ > 7^- Therefore, it might be interesting to consider also in the near future 
the possibility of having branches of the system under condition 7^ 7^ 7^- 
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(a) 




(b) 



Figure 2: (j)^ (a) and (jP' (b) as functions of ^ at the first four turning points of the spiral in 
the IV-characteristic. The lowest solution corresponds to the first turning point. The i-th 
solution is shifted {i — 1) unit(s) in the vertical direction for clarity. The space variable x 
has been scaled such that it is in the interval [0, 1]. 
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Figure 3: The IV-charactcristic of the stacked Josephson junction with J > 1. The critical 
current ratio J in this calculation is taken to be 1.5. Using AUTO, the 'back-bending' curve 
can be traced. It is known that solutions which lie on the upper/right branch are unstable. 
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Figure 4: Two pairs of solution to Q for J = 1.5 and 7 = 0.17 on the left and right 
branch of the IV-characteristic are compared. Picture (a) and (b) compare the 4>^'s and 0^'s 
respectively. The solutions on the right branch is given in dashed-dotted line. 
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Figure 5: The same as FigureOlbut with J < 1 where J is taken to be 0.5. In (a) we present 
the well known IV-characteristic where the curve can pass the lowest Swihart velocity c+. 
Because of the singularity, the AUTO calculation here stops at c+. In (b) another IV- 
characteristic is shown. All solutions lying on this 'back-bending' curve are unstable. 
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Figure 6: Two pairs of solution to |(3J) for J = 0.5 and 7 = 0.17 that lie on the left and right 
branch of the 'back-bending' curve (see Fig. IHJb)) are compared. Solutions on the right 
branch are given in dashed-dotted lines. 
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